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Studies of implementation of the simulation of the Bose-Einstein correlation effects 
in the JETSET particle generator are performed. Analysis of dependance of the 
one-dimensional correlation function parameters on the presumed boson source 
size reveals appearance of the effective new length scale, which limits applicability 
of the simple momentum-shifting mechanism, employed by JETSET to simulate 
the Bose-Einstein correlation. Two- and three-dimensional correlation functions 
are analysed as well and compared to the DELPHI data. 


1 Introduction 

Recent interest in profound studies of the Bose-Einstein correlations (BEC in 
what follows) was sparked by reports on possible influence of this ohenomenon 
on the measured value of the W boson mass in e+e" annihilationEstima¬ 
tions of the strength of this influence were done using the Monte Carlo particle 
generator JETSET0, which includes a phenomenologmal model for two-particle 
BEC simulation. JETSET (together with PYTHIAla) is so far the only e+e” 
annihilation event generator which accounts for BEC. It is known to reproduce 
well majority of experimental data, including some basic features connected 
to BEC. However, more sophisticated studies reveal not only certain discrep¬ 
ancies between emerimental and model results, but also self-inconsistencies in 
the model itself qH. Therefore, it is of particular importance to establish the 
extent of applicability of this model and to have a full understanding of its 
advantages and drawbacks. 

This work is devoted to studies of the built-in JETSET algorithm which 
is used to simulate BEC in hadronic decays of Z boson. The following section 
describes the physical process and corresponding observables. Overview of the 
studied model is presented in the next section. The last section contains results 
and discussion. 


2 Scope of the study and definitions 

The goal of this analysis is to study methods and consequences of implementing 
Bose-Einstein correlation models in the JETSET event generator, and to find 



out eventually how can it affect data analysis. To accomplish this goal, Bose- 
Einstein correlations between pions produced in hadronic decays of the Z boson 
are investigated. These decays are generated by the JETSET event generator, 
which uses the so-called Lund string model to simulate electron-positron anni¬ 
hilation events at a given center of mass energy, in our case Ecm = 91-2 GeV. 

Only the correlations between pairs of particles are studied, whith the 
two-particle correlation function is defined as 


C'2(P1,P2) 


P{Pl,P2) 
P{pi)P{p2) ’ 


( 1 ) 


where pi and P 2 are four-momenta of two particles, P(pi,P 2 ) ~ two-particle 
probability density, and P{pi) and P{p 2 ) denotes single-particle probability 
densities. 

Defined as above C 2 is parameterized in terms of the invariant four-momenta 
difference Q = ^ (pi - p-if - {Ei - £ 2 )^ as 

C2(Q) =iV(l + Ae-^''3=) . (2) 

This parameterization is one of the most commonly usedi, with the parameter 
R giving the width (or source size) and A - the strength of the correlation. 

In order to be able to compare our results with experimental data, only 
charged pions were used in the analysis. Although prompt pions produced in 
the string decay are the most clean sample, we selected only pions produced 
in decays of prompt p mesons, due to the following advantages : 

- high sample uniformity: no admixture of particles which are not subject 
to Bose-Einstein correlations; 

- sufficient average multiplicity of the sample: more pions are produced in 
decays of p mesons then promptly from the string. 

Studies were performed not only for the one-dimensional correlation func¬ 
tion (72(Q), but also for two- and three-dimensionaLcases. To facilitate cal¬ 
culations, the Longitudinal Centre-of-Mass System l3’B (LCMS) was used to 
measure four-momentum difference. This is the system in which the sum of 
the two particles momenta is perpendicular to the jet axis, hence only two-jet 
qq events were generated for further simplicity. In LCMS, Q is resolved into 
Qiongi parallel to the jet axis, Qt,out, collinear with the pair momentum sum, 
and complementary Qt,side, perpendicular to both Qiong and Qt,out- It is more 
convenient in some cases to use only two-dimensional picture with longitudinal 

(5|| = Qiong and perpendicular Q± = ^Qt^out + Qpside- Parameterization of 
C 2 for these two- and three-dimensional cases was chosen correspondingly as 

(72(Q±,Q||) = A^(l + Ae-^^^^-‘^ii^S) , (3) 





C2{Qt,OUt,Qt,side,Qlong) = N{1 + Xe Q^sideR^^ide QL„g-RLng). 

( 4 ) 

In high-energy physics experiments involving detectors, it is difficult to 
construct the product P{pi)P{p 2 ) from Eq.([|) due to the phase space limita¬ 
tions. Therefore it is often replaced by Po{pi,p 2 ), which is equal to P{pi)P{p 2 ) 
in a hypothetical case of absence of all the correlations. To make our results 
comparable with experiment, we must construct a reference sample correspond¬ 
ing to Pq{pi,P 2 )- Therefore, the measured two-particle correlation function is 
calculated as the double-ratio, using the event mixing technique 0 : 
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Here Ng^{Q) is number of like charged pions as a function of the four-momenta 
difference Q in presence of Bose-Einstein correlations. Subscript ‘‘‘BE, mix” 
denotes same quantity but with pairs of pions picked from different events. In¬ 
dices ^^noBE” and “noBE, mix” correspond to analogous quantities in absence 
of BEC (i.e., the simulation of BEC is not included into the event generation). 


3 Model description 

As it was already mentioned, JETSET is the only particle generator which al¬ 
lows and actually includes an algorithm emulating Bose-Einstein correlations. 
Recall that BEC is the quantum mechanical phenomenon, which has to appear 
during the fragmentation stage. However, in the standard implementation of 
BEC in JETSET, the fragmentation and decays of the short-lived particles like 
p are allowed to proceed independently of the Bose-Einstein effect. The BEC 
simulation algorithm is applied to the final state particles, for which the four- 
momenta difference Qij is being calculated for each pair of identical bosons 
i,j. A shifted smaller Q'^ is then to be found, such that the ratio C 2 {Q) of 
“shifted” to the original Q distribution is given by the requested parameteri¬ 
zation (Gaussian or exponential). In our case, the Gaussian parameterization 
identical to the form (||) was used : 

C2(Q) = l + (6) 

where A and Rinp are input parameters of the model. The input value of A is 
often set to 1, as it was done in this analysis too. Values of Rinp usually are 
chosen to fit experimental results. 

Further, under assumption of a spherical phase space, Q' is the solution 





of the equation : 
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After applying corresponding four-momentum shift to each pair of con¬ 
sidered bosons, all particle momenta are re-weighted to satisfy the energy- 
momentum conservation. This built-in JETSET algorithm works only in terms 
of the invariant four-momenta difference Q, i.e., it does not distinguish be¬ 
tween different Q components. This is yet another ambiguous assumption of 
the model. Also, it does not include particle correlations of higher orders. 

Evidently, this algorithm is absolutely phenomenological and is not based 
on any fundamental theory. It solely changes the final state particles momenta 
in order to resemble presence of the Bose-Einstein correlations. Moreover, 
presumption of the spherical shape for the phase space in Eq.(|^) is correct 
only for the case of very low Q (see the following discussion). 

In spite of all the ambiguities, JETSET reproduces fairly well experimen¬ 
tal data, such as shift of the p mass and observed Bose-Einstein correlations in 
terms of Q. It is widely used to calculate acceptance corrections for detectors 
and for various estimations, like the W mass shift mentioned in the Introduc¬ 
tion. To our mind, this peculiarity is worth investigating, if not in order to get 
better understanding of the BEC influence in experimental data, then at least 
in order to establish limits of applicability of such a simulation model. 


4 Analysis and discussion 

One of the most puzzling inconsistencies in the JETSET simulation of BEC 
is that the input shape of (|^) can not be obtained with the same parameters 
by fitting the resulting measumd <72 (Q) with formula (||) (see, for example, 
article by Fialkowski and WitB). This is mostly due to the improper phase 
space approximation in Eq.(^. However, this approximation can still be valid 
for certain input boson source size Rmp- To find out whether it is true, we 
studied 02 (Q) for different input values of Rinp in formula (^. 

Measured as the double-ratio (^ two-particle correlation function gener¬ 
ated with different input source size Rmp was fitted by the form (H). The A 
parameter always was reproduced at values close to 1, due to the high purity of 
the sample. The output source size R, however, behaved differently, see Fig. ^ 
and the corresponding table. Preliminary DELPHI results measured in 1991- 
1995 on all the charged particles are shown in the same table for comparison. 
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Figure 1: Width R of the measured in 

JETSET generated events correlation function 
C 2 {Q) as a function of the input source radii 
Rinp- Line represents the expected depen¬ 
dence. 
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Figure 2: Components of the correlation width as a function of the input source radii Rinp • 
a) i?|| , and b) Rt,out, ^tjSide^ ^tong • 
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0.459 ± 0.003 
0.616 ±0.005 
0.873 ±0.015 
1.530 ±0.070 

0.193 ± 0.001 
0.266 ± 0.002 
0.497 ± 0.005 
1.013 ± 0.021 
2.000 ±0.090 

0.562 ±0.005 
0.544 ± 0.004 
0.637 ±0.005 
0.814 ±0.011 
1.291 ±0.039 

0.499 ± 0.004 
0.511 ±0.004 
0.636 ± 0.005 
0.940 ±0.012 
1.622 ±0.051 

0.193 ± 0.002 
0.307 ± 0.003 
0.549 ± 0.004 
0.986 ± 0.013 
1.835 ± 0.065 

DELPHI 

0.300 ± 0.040 

0.640 ±0.020 

0.364 ± 0.009 

0.173 ±0.019 

0.628 ± 0.008 


It is clearly seen that the measured R does not depend on the input Rmp 
when the latter is below 0.6 fm. For the higher values of Rmp JETSET 
basically reproduces the demanded correlation function. 

Knowing that JETSET does not distinguish between components of in¬ 
variant momentum difference Q, we should expect similar behaviour of radius 
parameters of two- and three-dimensional correlation functions. Parameteriza¬ 
tion of these functions is performed in a form of a multi-dimensional Gaussians 
and (^) correspondingly. 

As one can see from Eig. ^ and in the corresponding table, transverse 
radii follow the same pattern as the R, while the longitudinal radius tends to 















reproduce the input value of Rinp- 

All these results show that there is a certain mechanism in the model, 
which imposes lower limit of around 0.6 fm onto measured R and its trans¬ 
verse components, and almost does not affect the longitudinal radius. The 
explanation of this phenomenon is illustrated at Fig. for the one-dimensional 
case. It shows evolution of the dN/dQ distribution with input source radius 
Rinp in comparison with the original non-correlated distribution. It is clearly 
seen that the expected Bose-Einstein enhancement appears only to the left of 
the non-correlated distribution peak, Q < 0.3 GeV. Since the model conserves 
multiplicity, and because the assumption of the spherical phase space in Eq. ( 0 ) 
is valid only for this region of the linearly increasing dN/dQ, a depletion ap¬ 
pears for Q > 0.3 GeV, which results in a non-Gaussian output correlation 
function. 



Figure 3: Distribution of Q as a function of the input source radii Rinp. Shaded area shows 
this distribution in absence of BEC. 


Therefore, the position of the peak in the non-correlated dN/dQ distribu¬ 
tion constitutes the limitation of the measured correlation width R and can be 
interpreted as a new length scale. This conclusion is also valid for transverse 
correlation radii (see Fig. In the longitudinal direction, dN/dQ^\ has less 
rapid falloff and peaks at a very small Q\\ value due to the LCMS properties, 
and thus is virtually insensitive to the mentioned length scale. 

As a result, one should state that the built-in JETSET model for simu¬ 
lating BEC is fully applicable only for sufficiently big sizes of boson source : 







Figure 4: Two-dimensional distribution d?N/dQ±^ dQ|| in absence of BEC. 

above 1 fm. At the same time, experimental data indicate that this size is 
around 0.5 fm at the Z peakla. This means that, at first, this kind of model 
has no predictive power. At second, one should be very careful when using 
JETSET with this model for calculation of detector corrections, because it 
can not produce an adequate unfolding matrix. 

It has to be mentioned that nowadays some other models for simulation of 
BEC are being developed 1131121. They use direct implementation of the Bose- 
Einstein interference into the string model, being theoretically accurate in this 
sense. At the moment they take too much computing resources to be used 
by high-energy physics experiments, but they do have predictive power. One 
of the most interesting predictions is that the transverse component of the 
boson source size, R±, must be significantly smaller then the longitudinal one, 
i?ll. As one can see, it is being confirmed by the preliminary DELPHI results, 
but it is not the case for the present JETSET version. This must encourage 
further works towards developing and implementing advanced models for the 
BEC simulation in particle generators. 
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